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We find a fine structure in the Aharonov-Bohm effect, wfiich is characterized 
by the appearence of a new type of periodic oscillations having smaller fractional 
period and an amplitude, which may compare with the amplitude of the conventional 
Aharonov-Bohm effect. Specifically, at low density or strong coupling on a Hubbard 
ring can coexist along with the conventional Aaronov-Bohm oscillations with the 
period equal to an integer, measured in units of the elementary flux quantum, two 
additional oscillations with periods 1/N and M/N . The integers N and M are the 
particles number and the number of down-spin particles, respectively. 

Prom a solution of the Bethe ansatz equations for N electrons located on a ring 
in a magnetic field we show that the fine structure is due to electron-electron and 
Zeeman interactions. Our results are valid in the dilute density limit and for an 
arbitrary value of the Hubbard repulsion U 

For the appearance or the existence of integer-fractional M/N flux quantum peri- 
odic oscillations the key role is played by the Zeeman energy. This is in an agrement 
with an assumption by Choy about importance of Zeeman energy. Without Zeeman 
energy there is a coexistence of the fractional 1/N and half-flux quantum periodic 
oscillations, only. With increasing magnetic field, the Zeeman energy increases and 
the half-flux quantum period transforms to the integer one. This transformation is 
not continious but via the appearance of fractional M/N oscillations, i.e. by some 
singular way. 

The fractional 1/N regime occurs when the value N < (LU/t)^^^, where L is the 
number of sites. The coexistence regime of integer-fractional M/N and fractional 
1/N AB oscillations occurs when {LU/tf/"^ « N « (LU/t). The parity effect 
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disappears in the fractional regime, but returns in the integer-fractional regime. We 
discuss the relation of the fine structure or the coexistence regime to existing exper- 
iments on single rings and on array of rings, where integer and half flux quantum 
periods are, respectively, observed. The explanation for the fractional 1/4 AB period, 
recently observed by Liu et al is proposed. 
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I. INTRODUCTION 



In the limit of large repulsion, when U/t ^ oo Kusmartsev |l| found that the Aharonov- 
Bohm(AB) effect is fractional. In units of elementary flux quantum its period is equal to 
where the is the number of electrons on the ring. This effect is drastically different 
from conventional integer or half-integer AB effects 0, 0, [Q, ||^ [^]. 

This result has since been confirmed in the investigation by Schofield, Wheatley and 
Xiang . The first corrections were studied by Yu and Fowler where they first recognized 
the importance of the magnon excitations for the Aharonov-Bohm effect. We extend their 
analysis and investigated the effects for the dilute system. We find a new type of oscillations, 
related to a magnon excitation spectrum. The new AB oscillations with the period M/N 
arising with the first corrections to the strong coupling limit or to the dilute density limit 
appears when we take into account the polarized or partially polarized ring with M down- 
spin electrons. That is for the effect it is very important a role of the Zeeman interaction 
causing the polatization effects. The importance of the Zeeman energy was noticed recently 
by Choy i. 

We show that for the very dilute system, when the filling factor tends to zero, as well 
as for strong-coupling limit, the fractional Aharonov-Bohm effect may coexist with 
the integer conventional, integer- fractional or half-flux quantum periodic AB effect. By the 
integer-fractional regime of AB effect we mean the appearance of a new type of oscillations 
on the envelope function of periodic oscillations, which display a new fractional period. 
The period of these oscillations is equal to M/N, where M is the number of down-spin 
particles. The M/N oscillations exist even when the oscillations are washed out. When 
the Zeeman interaction is neglected, i.e. the Zeeman energy is zero and M = N/2, there is 
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a coexistence of the fractional and half-flux quantum periodic oscillations. In contrast, 
when the system is almost polarized, i.e. M = 1, the fractional and integer periodic 
oscillations coexist. 



In previous work |T0| we found that finite scaling finite size effects were very important. 
There is a scaling behavior of the ground state energy which does not depend on the size L 
or on f/, but depends only on N/UL = const = a, where U is measured in units of t and L 
is measured in units of the lattice constant. Such scaling occurs only at small values of a. 
Due to this scaling symmetry the fractional Aharonov-Bohm effect may arise even for small 
values of U, for the very dilute electron systems. In the present work, we shall show that, 
such a symmetry is related to the structure of the low energy excitations of the Hubbard 
Hamiltonian, namely, with the spectrum of the spinon or magnon excitations. In fact the 
parameter a defines the energy of the magnon excitations. Thus, the shape of the new AB 
oscillations completely depends upon the energy of the magnon excitations. 

In fact, we get an amusing physical picture. First, let us discuss the case M = N/2. For 
an odd number of particles in the ground state of the Hubbard ring, there already exists 
a spinon excitation. The momentum of this excitation is associated with the top of the 
spin- wave spectrum. This excitation is needed because of the parity effect, discussed by 
Legget Loss and Kusmartsev for the case of interacting spinless fermions. 

Thus, again with the first correction to the strong coupling limit or to the dilute density 
limit we have shown that on the Hubard ring there occurs the reconstruction of the parity 
effect, reflecting the structure of low-energy excitations of the Hubbard ring. However, the 
parity effect obtained is distinct from the parity effect for the mesoscopic ring with spinful 



electrons and disorder |T3| Inducing the flux the momentum of this spinon excitation, which 
exists in the ground state, changes. This gives rise to the correspondence between the spin 
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excitation spectrum and the Aharonov-Bohm fine structure. With the magnetic field the 
number of polarized electrons changes, i.e. the spin excitation spectrum changes, this gives 
rise to the change of the Aharonov-Bohm periodicity. 

This occurs also for the general situation, with any number of particles on the ring and 
at different rational ratio M/N . In the latter case one may see an additional fractional 
M/N AB oscillations, which are due to the structure of the spinon excitation spectrum; the 
picture is as follows. With the fiux, the spinon changes its momentum, which is, in fact, 
proportional to the fiux. The change in the momentum does not occur continuously, but via 
jumps related to the finite number N of the particles on the ring. Therefore these jumps 
are associated with the fraction of elementary fiux quantum. Thus, on the Hubbard 
ring with the fiux the spinous behave similar to spinless fermions, where the momentum is 
also shifted with the fiux, reconstructing the parity effects. The Fermi momentum of this 
fermions is related to the number of polarized electrons M. This gives the new scale and the 
new periodicity in Aharonov-Bohm effect, related to the ratio M/N. However, the fractional 
M/N AB effects do not destroy the fractional effect associated with the bound state of 
N particles. Both effects may coexist, refiecting the separation of spin and charge degrees 
of freedom. 

II. MAIN EQUATIONS 

To demonstrate this picture we have studied the Hubbard Hamiltonian 

L 

H^-t al^aj^ + Uj2iT-i+iT-i- (2-1) 

<i,j>,<T i=l 

which is parametrized by the electron hopping integral t, the on-site repulsive Coulomb 
potential U and the number of sites L. The operator al^{aia) creates (destroys) an electron 

6 



with spin projection a {a = +or— ), at the ring site 2, and njo- is the occupation number 
operator a\^ai„. The summation in Eq. (1) extends over the ring sites z or - as indicated by 
< i,j >,a - over all distinct pairs of nearest-neighbor sites, around the ring with the spin 
projection a. For the case of the magnetic field we use the same form of the wave function 
as in Refs 0, and [!§], i.e. 



N 



J2[Q,P] exp[iJ2kpjXQj] 
p j=i 



(2.2) 



where P = (Pi, . . . , P/v) and Q = (Qi, . . . ,Qn) are two permutations of (1, 2, A^), and 
N is the number of electrons. 

The coefficients [Q,P], as well as (/ci, .../ctv), are determined from the Bethe equations, 
which in a magnetic field are changed by the addition of the fiux phase 27r/ ||16|, [|l^], [jl9 



1^, 20 , 21 , and 22 



/3=i ^ sin kj -iXp + U /4 



(2.3) 



and 



^ it sin kj — iXa — U/A 



^ - iXp + U/2 
^},\^K-^X^-U/2 



(2.4) 



fJl^it sin kj — iXa + U/A, 
f being the fiux in units of the elementary quantum fiux 0o- The explicit form of the Bethe 
equations in a magnetic field is [[1^, [|^, [|| 



M 

Lkj = 2nlj + 27r/ - ^ ^(4(t sin kj - A^)/f/) 

/3=1 

N AI 

Y.e{A{tsinkj- Xp)/U) =2TxJp+ 0i2{Xp-X^)/U) 



(2.5) 
(2.6) 



where 6{x) = 2arctan(a;) and the quantum numbers Ij and Jp, which are associated with 
the charge and spin degrees of freedom, respectively, are either integers or hal- odd integers, 
depending on the parities of the numbers of down and up-spin electrons, respectively, i.e. 



M 

Ij = — (mod 1) and Jp 



N-M+1 



[mod 1). 



(2.7) 



III. EXPANSION 



In previous work showned that in the dilute density hmit there exists a small pa- 
rameter a, which is 

« = ^ = Pt/U, (3.8) 

where p = N/ L is the filling factor. 

From an expansion in a ( a << 1) of the second Bethe equation (|2.6|), we get the 
equations 

8 1^ Af 

^^^^^^ - 777TT72T E Mkj) = 2nJ^ + 0{{tp - t,)/2), (3.9) 

U -t- j=i A„=l 

where t/j = AXp/U . 

Next step, we notice that the second term on the left hand side the eq.( |3.9| ) is small and 
therefore we may incorporate it into the argument x of 9{x) — , giving 

iV^[t^-^Esin(fc.)J =2vrJ/5+ E e{{tp - t^) /2), (3.10) 

With the substitution xp = tp — I^jLi sin(A;j), this equation is reduced to the Bethe 
equation for an isotropic Heisenberg antiferromagnet on a ring having sites and M down 
spins: 

M 

N9{xp) = 2nJp + e{{xp - x,)/2), (3.11) 

a=l 

One sees that the solution for Xa is independent of both the fiux / and on its repulsion U. 
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Let us substitute the variables Xfs = U{xp + SjLi sin(fcj))/4 into equation (|2.5| ) for 
the momenta kj and study the resulting equation in the variables xp. Then, analogously to 
the expansion of eq. (^^ ), expanding the eq.( |2.5|) in powers of a , we get a new form of the 
first Bethe equation, namely 

f3=l /3=l-L+^/3 

For parameter a small (dilute density), this system of linear equations may be solved, with 
the result for the momentum 

n -I M RR 1 ^ -I M 

and for the ground state energy 

E,rouna = -D C0S(| (/ - ^ + ^ + ^ E + / - ^ )) ) (S-^) 

where p = — Y^^Li Jp , Imax and Imin are the maximal and minimal charge quantum numbers, 

b = 2sin(7rA^/L)/sin(7r/L), 

and 



The magnitude 



M ^ 



/3=1 -^p 



which is a real number, has the physical meaning of the total energy( ground state + 
excitations) of the antiferromagnetic Heisenberg Hamiltonian with the exchange constant 
equal to 1/4. 
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The variables Xa are satisfied by the Bethe ansatz equation for that Hamiltonian ( P . 1 ID ■ 



Actually, the value B{k) is a function associated with the excitations spectrum, where k is 
the momentum of the spin-wave excitation of the Heisenberg Hamiltonian k = 27Tj2a Ja/N. 
For the antiferromagnetic state of Heisenberg Hamiltonian the dependence B{k) has been 
studied in detail in many works: first being estimated by Anderson , within a spin wave 



approximation and then by des Cloizeaux and Pearson in the framework of the 
Bethe ansatz obtaining the exact form of this spectrum. 
At large number N this dependence has the form: 

B{k) = {N{\og 2 - 1/4) - ^ I sin A; |)/4 (3.15) 



The ground state energy(first two terms) has been obtained by Hulthen p6|] in the framework 
of the Bethe method It was seven years later after Bethe introduced his method (the 
Bethe ansatz), applying it to the ferromagnetic Heisenberg chain. 

IV. SINGLE SPIN 

Let us consider the case with M = 1, i.e. when the system is strongly polarized. Here, 
we find an analytical solution, permitting eq. (|3.11| ) to be immediately solved yielding 



x, = tan(^) (4.16) 

The function B{Jfs) can then be calculated from this solution, giving B = cos^(^^). 

Now we must consider two independent cases, when is odd and when A^ is even. Let 
A^ be an even number, then is an integer, p say. Then, for the ground state energy we 
get the explicit expression 
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The ground state energy-flux dependence is this a function with two types of oscillation. 
The period of the first of oscillation is equal to the single flux quantum, while the period of 
the second type of oscillations is equal to the fraction 1/N oi the elementary flux quantum. 
Let us describe this dependence. Each period of the single flux quantum then consists of 
A'"-parabolic like curves, where each of these curve is labelled by p: < p < N. In the hmit 
of large N we may put as f = p/N to get an exact single flux periodic function, which has 
minima and maxima at integer flux values and at half-odd integer flux values respectively. 

For an illustration on Fig.l is presented the energy-flux dependence for the Hubbard ring 
having 1000 sites with 6 electrons and a single down spin at the coupling U/t = 10. On this 
Figure one sees 7 parabolic curves having minima at p/6 flux values, and 6 cusps associated 
with intersections between each two neighbour parabolas. The envelope of these parabolic 
curves is the single flux periodic function cos^ tt/. 

Next we turn to the case of N odd, when the numbers Jg are half-odd integer. So, let 
us take Jg = p + N/2, where p is an arbitrary integer. Given this, then the ground state 
energy may be expressed as 

^ L+-^ sin''(j>n/N) ' 

Again, we find that the ground state energy dependence has both the single fiux quantum 
and the fractional flux quantum periodicity, described by the flrst and the second factors, 
respectively. In the region of the one period of the single flux quantum this dependence 
again consists of iV-parabolic like curves, labelled by the integer p: < p < N. In the limit 
of large N we may put as f — p/N to get an exact single flux periodic function. However, 
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for an odd number of particles, this function now has maxima at integer flux values and 
minima at half-odd integer flux values. 

For an illustration on Fig. 2 is presented the energy-flux dependence for the ring having 
1000 sites with 11 electrons and a single down spin at the coupling U /t = 10. On this Figure 
one sees 12 parabolic curves having minima at p/11 flux values, and 11 cusps associated with 
intersections between each two neighbour parabolas. The envelope of these parabolic curves 
is the single flux periodic function ~ — sin^ vr/. 

Several conclusions can be drawn immediately. Firstly, the parity effect, which disap- 
peared in the a small limit, is again recovered when the first correction is taken into account. 
In other words, the ground state energy-flux dependence is different for the even and odd 
numbers of fermions on the ring. To be more precise, the behavior of the flux-energy de- 
pendence for an even number of particles is shifted by a half flux quantum in comparison 
with the case of an odd number of particles. This is similar to the parity effect observed on 
a ring with interacting spinless fermions |Tl]|, |T^. However, in comparison with spinless 
fermions on the Hubbard ring the, cases of even and odd number of particles are reversed. 
The parity effect discussed might be related to the number of down-spin electrons, which 
is — 1. In comparison with the case of spinless fermions, because of this shifting of the 
particle number by — 1, the parity effect discussed might be reversed. In this case, 
again, one may connect the parity effect with the creation of the statistical flux vr, due to 
Fermi statistics [^]. 
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V. SPIN UP-DOWN EQUILIBRIUM 



In this section we will follow the important suggestion by Yu and Fowler about the 
possible role of spin excitations in the Aharonov-Bohm effect on the Hubbard ring in the 
limit of large U/t. Our results extend their analysis of the first corrections in the strong 
coupling limit as well as being applicable to systems with the low electron density at arbitrary 
coupling. Now, let us consider the general case, when the number of spin-up fermions is 
approximately equal to the number of spin-down fermions, i.e. it may differ by ±1. 

To investigate the behavior of the fiux-energy dependence for such a general situation, 
we must consider four independent cases, when the numbers and M are even or odd. As 
we will see below ( also, see, above the single spin case M=l) for the ground state energy 
flux dependence only the parity of the number A^ will be important. 

Let us consider the first of these cases, when A^ and M are even. The quantum numbers 
(which are half-odd integer) of the Heisenberg Hamiltonian, J^, for the ground state are 
distributed symmetrically with respect to origin, so that Jp = 0. The quantum numbers 
of the charge degrees of freedom 1^ for the Hubbard Hamiltonian are integer with J2 = 
-f = -M. 

However, at zero flux the minimum energy of the Hubbard Hamiltonian will correspond 
to the set of quantum number { Jq}, but with nonvanishing sum Y.Ja = M, which does not 
correspond to the ground state of the Heisenberg Hamiltonian. Such a state is associated 
with the single magnon excitation of the Heisenberg chain, having a momentum ko = = 
vr and the excitation energy B{ko) = Bq. The distribution of quantum numbers 
associated with the ground state of the Hubbard ring at zero flux has a hole near the origin, 
i.e. is of the form: 
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With increasing flux this hole in the distribution is shifted to another place so that Ja = 
M — 1 and, therefore, the magnon momentum decreases. In other words with the flux the 
hole in the distribution moves to its left side. Therefore, with that motion of the hole, 
exactly speaking, with the new position of the hole the magnon excitation will have the new 
momentum k < kQ and the energy B{k) < Bq. Exactly speaking, the obtained dependence 
B{k) describes the low-energy excitation of the discrete antiferromagnetic Heisenberg chain. 



as first discussed by des Cloiseaux and Pearson |Q p5[| . 

With the dependence B{k) taken into account the appropriate ground state energy de- 
pendence of the Hubbard model on the flux may be written in the form: 

Et:::z^' = -d cos [^(/ - ^)], (5.20) 

where the flux value / is changed in the region (2p — 1)/2N < f < (2p+ 1)/N and D{k) = 
2 sin{TT N/L{k))/ sm{TT/L{k)), with k = 2Tip/N being the magnon momentum of the finite 
Heisenberg chain with sites, where p takes integer values. 

The dependence consists of equidistant parabolic segments. The position of each 
parabolic segment along the vertical axis depends upon the value —D{k). As the inte- 
ger p increases, the magnon momentum k and the energy of this magnon excitation B[k) 
also increases, so long as p < y. The magnon energy then has a maximal value at p = y- 

This means that the value —D{k) associated with the position of bottom of the p-th 
parabolic segment along the vertical axis decreases monotonically. As a result, we see that 
the envelope curve which passes through the bottom of these parabolic segments is similar 
to the curve describing the spectrum of the spin-wave excitations of the antiferromagnetic 
Heisenberg chain with A^ sites. Note, that this curve is described by the equation Egr{k) = 
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—D{B{k)). The flux dependence of tliis envelope curve can be obtained if we substitute 
k = 27r/. It is a half flux quantum periodic function. 

Now, let us consider the second case when is even, but M is an odd number. In this 
case, at zero flux, the ground state energy corresponds to half-integer numbers Ij and integer 
numbers Ja distributed symmetrically around the origin (see, for example 0). That is, here, 

E^.=E'^/3 = (5.21) 

The analysis of this case is similar to that above. The given distribution of quantum 
numbers Ja corresponds to the ground state energy, both of the antiferromagnetic Heisen- 
berg chain and the Hubbard ring. However, including the flux the ground state energy of 
the Hubbard ring will now correspond to the nonvanishing sum J^/s Jf3, which is negative. 
This means that there will be created a magnon excitation similar to that described above, 
associated with the appearance of a hole in the distribution of Ja- With the increasing flux 
the hole will move from the left to the right of the distribution. 

We, therefore, once again have the flux-energy dependence of the Hubbard ring de- 
termined via the low energy excitations of the Heisenberg chain. This dependence again 
consists of parabolic segments and is described by a similar formula to eq.( |5.20D . The 



envelope, which passes through the bottom of these parabolic segments is described with 
the aid of the same equations —D{27[f), i.e. it has cusps at its minima, which are at integer 
and half-integer flux values and smooth maxima at the values of the flux exactly midway 
between the cusps. 

Let us now consider the case when is odd and M is even. The quantum numbers Ij are 
integer and are symmetrically distributed around the origin, i.e. = 0. The quantum 

numbers Ja are also integers. The ground state energy of the Heisenberg Hamiltonian 
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corresponds to the dense distribution with the sum Y^aJa — ~~2- However, the ground 
state energy of the Hubbard ring with zero flux values must correspond to the state with 
J2a Ja = 0. This sum can equal zero when the associated state of the Heisenberg chain is 
not a ground state, at rather an excited state. This means that it is a state with the hole 
in the {J a} distribution. One sees that in this case the hole will be created exactly in the 
origin: 

M ^ ^ M , , 

Ji, . . . , Jm := - y, . . . , -1, hole, 1, ■ ■ ■ , y (5.22) 

This state corresponds to the top in the spin-wave excitations of the Heisenberg chain, 
i.e. to the spin- wave excitation with momentum k — ^. As the flux changes, the parabolic 
segment with its minimum at zero flux is exchanged with another the parabolic segment 
which has its minimum at the flux f = jf- 

This parabolic segment will correspond to the new set of J^, where the hole is in the 
position formerly occupied by = 1, so now J2a J a — ~1- The appropriate magnon 
excitation associated with this parabolic segment will have momentum k = and 
energy B C^^^'^^ ), which are lower than the momentum k — tt/2 and the maximal energy 
S(7r/2), respectively. 

With the next increasing the flux, there will occur a transition to a new parabolic seg- 
ment, corresponding to a new set of Jq, with the hole at the position — 2. This is the 
parabolic segment corresponding to the magnon excitations with momentum k — 
which has a lower energy than the magnon excitation in the previous case. Therefore, at 
this flux the mimimum energy of the appropriate parabolic segment decreases. The absolute 
minimum is reached at the quater and 3/4 of the flux quantum, if the flux is changed within 
a single flux quantum. 
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The envelope curve, which goes through the minima of these parabohc segments is a half 
flux quantum periodic function, which has smooth maxima at integer and half integer flux 
values. This curve may be described with the aid of the formula E ground = ~-D(27r/). 

Let us consider now the last case, which is N and M odd. Here, both the and Ij are 
half-odd integer. The energy of the Hubbard Hamiltonian, which corresponds to the dense 
distributions with Y^Ij — ~y and with J2Ja=^ and assumed to be a ground state is not 
a ground state of the Hubbard ring at zero flux. But this dense distribution of { Ja} is the 
true ground state of the Heisenberg Hamiltonian. 

The ground state of the Hubbard ring will correspond to the distribution of the quantum 
number {Ja} with the hole at the position — and the sum Y.^=i — +^^^|^ ~ +T- 
That is, to say, in the ground state of the Hubbard ring there already exists the magnon 
excitation with momentum k = ^^^^tt = 7r/2, which is, in fact, associated with the top 
in low energy spin wave excitation spectrum of the Heisenberg chain. The energy-flux 
dependence will again consist of parabolic segments as described above. The minima of 
these segments will be located at the flux f = ff- The energy behavior of the minima of 
these segments will be exactly the same as described in the previous case. 

The picture discussed is illustrated in Figs 3 and 4. In Fig. 3 one sees the coexistence 
of the half- integer and of the fractional 1/A^ Aharonov-Bohm effects. The ring has 7500 
sites filled with 25 electrons. There are two types of oscillations present: half-flux quantum 
periodic, which has the largest amplitude and the other fractional - 1/25— flux quantum 
periodic, the amplitude of which is much smaller. The envelope curve of the fractional 
oscillations is a half flux quantum periodic function, with smooth maxima at integer and 
half-integer flux quanta / = —1/2, 0, 1/2, 1, ... and cusps, which are minima, in between. 
This is due to the parity effect, since the number of particles on the ring is odd. 
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For contrast, in Fig. 4 we present a case for an even number of electrons on the ring with 
= 10, and L = 1000. Again, as in Fig. 3 the coexistence of the half-integer and fractional 
Aharonov-Bohm effects is observed and because there are 10 particles the fractional 1/10 
flux quantum oscillations arise. The envelope curve of these fractional oscillations is as in 
Fig. 3 a half flux quantum periodic function, which, however, has minima at integer and 
half-integer flux quanta. 

In the comparing these two Figures, one notes that with the change of the number of 
particles from odd to even, the position of the maxima is shifted by a quarter-fiux quantum. 
This is precisly due to the parity effect similar to the parity effect for spinless fermions 
0, |jl2[, [|lT|. That is on the Hubbard ring we have obtained the parity effect which is 
similar to one for spinless fermions. The effect is due to the contribution of the spinon 
excitation spectrum in the AB effect. Note, that in the second case (Fig. 4) the difference 
in amplitudes between the half-integer and fractional 1/A^ oscillations is smaller, indicating 
that the fractional effect is more important at small number of particles. With increasing 
U or L, this difference may disappear. 

VI. FRACTIONAL ^ AHARONOV-BOHM EFFECT. 

In previous sections we have shown that the AB periodicity for the cases of M = 1 and 
M = N/2 are different ( single and half flux quantum, recpectively). Therefore, there occurs 
the question: what kind of AB periodicity occurs at arbitrary M: 1 < M < N/2 Let us 
now consider the case when M is arbitrary, but M « N . For this general case, when the 
density of down-spin particles M/N is much smaller, then the density of up-spin particles 
1 — M/N can again be found an analytic solution. In this case we expand equation (p.ll|) 
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with in the parameter M/N . We can assume that the values of xp in eq. (|3.11|) are smalL 
Using these assumptions the eq. may be written in the form 

X, + V ^""^^ = ^ (6.23) 

i.e. we have obtained a system of hnear equations, which can be immediately solved, giving 

_ _ ^ irJq 

^^~iV-2M ^^{2N-M)N ^ ' 

This solution looks like the momentum spectrum of a free particle on a chain of 2N — M 
sites ( the first term of the right hand side of the equation), plus an additional fictitious flux 
associated with the second term of the right hand side of the equation. From the solution 
obtained one sees that the values xp are small if tiM/(2N — M) << 1, which is equivalent 
to our small parameter M/N << 1, and therefore << 1 does, indeed, hold. 

On the other hand this analysis has shown that the second term of the right hand side 
of eq. (p.ll ) is smaller that the first one in the parameter ^ <C 1. Therefore, in that case 
the solution may be obtained in the form: 

= tan (6.25) 

The substitution of which into -B[xq,] gives 

Bik) = - + -J2^cos^ (6.26) 

The choice of the set Ja determines the dependence of the low lying excitations B{k). Now 
we must consider four independent cases associated with the parity of numbers and M. 

Let, first us consider the case, when N is even and M is odd. Then the ground state 
energy corresponds to the vanishing sums of the integers Ja, J2a Ja = ^ and of the half- 
integers = 0. Both of these sets are densly distributed, symmetrically around the 
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origin. The excitation will correspond to the hole in the distribution of J a- For example, 
when p = 1 = —YjJai then the hole is created on the left hand side, from the left edge . 
For the value p = 2, the hole is moved on to the next position . In general, if the hole in the 
position J a then we can write 

M+l 

P=-T.Jp + J^ (6-27) 

/9=1 

where the sum is taken over the dense distribution. Therefore, if p < M, this sum equals 
YJp=i Jf3 = ~{^^^2^)- In the similar way, by splitting into to terms, the sum over the dense 
distribution of M + 1 quantum numbers J/j and the term associated with the hole at the 
place of Ja we find for the energy of such excitations: 

Then, expressing from eq.( |6.27| ) and substituting it into eq. ( |6.26| ) we get the dependence 
of B(^) = B{k) where p < M . Note, that the case of p = — S^li Jp = M corresponds 
to the ground state distribution of Ja shifted by 1 to its left side. When p > M, in this 
shifted distribution there will again appear the hole, initially, at the beginning on the left 
hand side, from the left edge. 

With increasing p this hole will move from the left to the right side, once M < p < 2M. 
At the value p = M, the distribution of will be dense, but shifted by 2 to the left in 
comparison with the distribution associated with the value p = 0. Thus, for an arbitrary 
value of p, we can write 

p = +liM + I) + {^^^) - Ja (6.29) 
The energy of the associated excitations is equal to: 
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r./ X Ir,. sin^^^^^ ,27r/ + 7r, 27r , M + 1 

B{p) = ^[M + _-Lcos(^^) - co5-(p - ^— - 1{M + 1))] (6.30) 

where / = 0, if < p < M, / = 1, if M < p < 2M, / = 2, if 2M < p < 3M and so on. 
Thus, the ground state energy has cuspoidal minima (cusps) at p = nM, where n is an 
integer. So, if we substitute p = IM, we have for the envelope functions of these cusps: 

Since p < y, the function S(Z) is monotonically decreasing with having a maximal value 
when / ~ [2^]) where the brackets indicate that we take the integer part of the value The 
spectrum of the excitations of the Heisenberg chain reflects the flux quantum periodicity 
of the Hubbard ring, as described above. Therefore, we come to the conclusion, that in 
addition to fractional 1/N, and integer on the Hubbard ring there also occurs fractional ^ 
flux quantum periodicity. 

For an illustration we give several examples for the ring with U/t — 10 and L — 10^, 
presented in Fig. 5. In the first example (see, Fig.Sa) the energy-fiux dependence is given for 
N — 20 and M = 3 in the region of the half flux quantum. Since the number of particles is 
not large one sees, the very pronounced 1/N oscillations associated with the single parabolic 
like segment. However, one sees already M/N oscillations with much smaller amplitude. So, 
the first, fourth and seventh parabolic like segments have minima lower than nighbouring 
parabolas. If the number of particles increases, for example, = 30 (see. Fig. 5b, where we 
taken M = 5) the amplitude of M/N oscillations increases. In this case the amplitudes of 
1/N and are M/N oscillations are, already, compared. With the next increasing number of 
particles N, for example, when = 40 and M — 7 presented in Fig. 5c, the amplitude of 
M/N oscillation becomes larger than the amplitude of the 1/N oscillation. 
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Let us consider the second case, when N is odd and M is even, which corresponds to two 
types of integer quantum numbers Jq, and Ij, with Z^/j = and J2a J a = However, 
the ground state of the Hubbard ring will correspond to an excited state of the Heisenberg 
chain with the hole in the origin, i.e. the ground state corresponds to the sum: J2a Ja — 0- 

The excitation associated with the value p—1 corresponds to the shift of this hole to the 
position Jq, = 1. So with the increasing p the hole moves to right. Therefore, we can write 
p—— Y,a=i Ja + J/s — Jf3: whcrc p < M/2 and the sum is taken over the dense distribution. 

As p is increased from M/2 the hole will be created on the left edge and moved to the 
right edge, where p < (M + M/2). With the next increase of p this process will be repeated 
and so on. Thus, we may write in general form that 

p = +l(M + 1) + Jfs, (6.32) 

where the first term related to the sum J2^=i Ja- 

The envelope function of parabolic like curves associated with the ground state energy- 
flux dependence is calculated in a similar way and equals: 

B{p) = -[M + ^^^1 cos{ — ) - cos-{p - 1{M + 1))]. (6.33) 

where / = 0, if -(M)/2 < p < (M)/2; / = 1, if (M)/2 < p < (M)/2 + M, Z = 2, if 
(M)/2 + M < p < (M)/2 + 2M and so on. 

This dependence is shifted by ^ is comparison with the previous case. It consists of 
cuspoidal minima and smooth maxima in midway between. In comparison with the previous 
case here the cuspoidal minima, one of which existed at the value p — 0, are shifted by 
the same value of the flux, too, and are exchanged in the positions with smooth maxima. 
Thus, the properties of the envelope function are the same as in the previous case. For an 
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illustration in Fig.6 we present the energy- flux dependence for the ring with U/t — 10. So, 
it is clearly coexist 6/19 ~ 1/3 and 1/19 flux quantum periodicities for the ring with 19 
electrons having 6 down-spins for number of sites L = 1000 (see, Fig.6a). The parabolic like 
curves are associated with 1/N type oscillation. The envelope function of these curves is 
6/19 ~ 1/3 flux quantum periodic function. The cusps of this function are located at 3/19 
and 9/19 flux quantum. With the increasing number of particles the amplitudes of both type 
of fractional oscillations decrease. This can be seen in Fig.6b, where the particles number 
was changed to 49 and the number of sites was taken as 10^. Although with the next increase 
of the number of particles at the flxed value of M the amplitudes of both type of oscillation 
decrease, the fractional 1/A^ oscillations begin to wash out, however the fractional M/N 
oscillations are still exist. This is shown in Fig. 6c where the number of particles was taken 
as = 99 on the ring with the same the number of sites, as in Fig. 6b. On the other hand 
when both numbers, the number of particles N and the number of dowm-spins M increase, 
the amplitude of M/N oscillation increases and it is compared with the amplitude of the 
main single flux quantum periodic AB oscillation. In the latter the fractional 1/A^ oscillation 
is completely washed out, although it is an origin for the fractional M/N oscillation. The 
fractional 23/99 AB oscillation is presented in Fig.6d, where we consisder the ring with 
M = 23 and A^ = 99. One sees on this Figure that the fractional 1/99 oscillation is 
practically washed out. 

The other two cases only differ from these already considered ones only by an appropriate 
shift. In both of them the ground state of the Hubbard ring corresponds to the excited state 
of the Heisenberg chain with the value 

M 

a=l ^ 
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Because of the M-periodicity, this state may be mapped into another one, where the 
sum '}2^=i Ja = Q vanishes, i.e. this can be obtained with the aid of the transformation: 

Ja Ja ^1) 

where h is the integer part of the ratio j^. Therefore, we can write 

a a 

where q < [y]- 

Thus, the problem has been reduced to the one, discussed above, but with the differ- 
ence, that the ground state of the Hubbard ring will correspond to the excited state of the 
Heisenberg chain where the hole in the dense distribution of the set Ja is located at the 
— g-th position, which is determined by eg. (|V^) . The other excited state of the Heisenberg 
chain (or the ground state energy-flux dependence of the Hubbard ring) is associated with 
the shift of this hole to the right edge. Then all excitations may be described exactly by the 
same manner as above. That is we may write 

M+l 

P=- E Jc. + Jp (6.34) 

a=l,dense 

Now we must consider the two particular cases ( 1) and M are even and 2) and M 
are odd), independently. 

In the flrst case, the values Ja or Ja are half-odd integer and from eq. (|6.34| ) we get 



M+l ~ , , 

P = + Jp- (6.35) 

The excitation energy is estimated as above yielding: 
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From eq.( |6.35|) and making the inverse transformation from J a to J a we finally get 



where the value / is choosen thus / = 0, if -N/2 <p< -N/2 + M; / = 1, if -N/2 + M < 
P < -N/2 + 2M, I = 2, if -N/2 + 2M < p < -N/2 + 3M and so on. 

For an illustration we present in Fig. 7 the dependence of the ground state energy for the 
ring with 1000 sites and U/t = 10. The number of electrons is taken as = 20 amd M = 6. 
One sees on this Figure that the dependence of the envelope function of the parabolic curves 
is monotonically decreased with the flux increases from zero. This dependence has cusps, 
which are repeated with the periodicity 6/20. So we have here a quasi ~ 1/3 flux quantum 
periodicity. 

Similarly way we can obtain the excitation energy for the second case, when the parity 
of particles number is odd, i.e. N and M are odd: 

B{p) = \[M+ -^—^cos'-^ - cos'-^iP + ^-m + 1)]] (6.38) 

where the value / is chosen thus / = 0, if -A^/2 - {M -l)/2 < p < -N/2 + (M - l)/2; 
/ = 1, if -A^/2 + (M - l)/2 < p < -A^/2 + (3M - l)/2, / = 2, if -A^/2 + (3M - l)/2 <p< 
-N/2 + (5M - l)/2 and so on. 

For an illustration we present in Fig. 8 the dependence of the ground state energy for the 
ring with 1000 sites and U/t = 10. The number of electrons is taken as = 19 and M = 5. 
One sees on this Figure that the dependence of the envelope function of the parabolic curves 
is monotonically decreased with the flux increases from zero. This dependence has cusps, 
which are repeated with the periodicity 5/19. So, here again as in the previous case we have 
a quasi- ~ 1/3- flux quantum periodicity. 
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Thus, for all four cases the excitation energy is as M/N periodic function, having cusps 
repeated with the same periodicity. To study the flux dependence we must substitute the flux 
value / instead of the value p/N oi the excitations of Heisenberg ring. Generally speaking, 
in the last two cases the ground state energy-flux dependence of the Hubbard ring related 
to these excitations of the Heisenberg chain is shifted by half-flux quantum in comparison 
with two previous cases, respectively. But, the qualitative shape of these dependencies are 
not changed. 

VII. FINE STRUCTURE 

One general conclusion that can be drawn here that the complete picture of the 
Aharonov-Bohm effect must include a fine structure. The latter appears when we take 
into account the electron-electron interactions and is characterized by oscillations with an 
amplitude and a period smaller than the conventional one. In fact, there is a hierarchy of 
flux periodicities which depends implicitly on the Zeeman interaction. 

The Fourier spectrum of the energy-flux dependence for the Hubbard ring mainly con- 
sists of two or three harmonics. The number and type of these harmonics depend on the 
contribution from the Zeeman interaction. One of these harmonics is conventional and has 
a single flux quantum period; the second one has the fractional 1/N flux quantum period; 
the third one has half-flux quantum period or fractional M/N period. The latter appears 
when the number of down spin particles M and the number of up-spin particles N — M are 
different. 

Let us analyse the fine structire. The amplitude of the single-flux or the fractional 
M/N periodic oscillations is proportional to the parameter n^cc, where n — 'kN /L and we 
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use dimensionless units. The amplitude of the fractional M/N oscillation is smaller, by 
factor M/N, than the amplitude of oscillations with the integer period. Here, we assume 
that M << N. When M ~ N/2 the situation is reversed and the M/N oscillation has a 
dominant amplitude. The amplitude of the fractional 1/N flux quantum periodic oscillations 
is proportional to ^ and does not depend on the value U. One sees also that both types 
of oscillations only depend upon U and L through their product [/L or in other words they 
only depend upon a. With increasing N at constant n and a, the amplitude of the fractional 
1/N oscillations vanishes. 

The single flux quantum harmonic coexisted with fractional M/N periodic harmonics 
show the parity dependence. That is there is a difference in the energy-flux dependencies for 
the cases of an even and an odd number of fermions. This is in contrast with the fractional 
periodicity, which has no the parity dependence. In the present case the parity effect 
for an even and an odd number of fermions is transposed with respect to the case of spinless 
fermions. 

It is useful explicitly to write an envelope of the flux- dependence. For example in the 
case, where there is up spin- down spin equilibrium (i.e. M=N/2), for the two cases of an 
even or an odd number of particles on the ring the flux energy dependence has the form: 



For the other values of M the pref actor will be the same, but the oscillation factor will 
be changed to reflect the fractional M/N flux quantum periodicity. 

On the other hand the amphtude of the fractional 1/N oscillations is estimated to be 
~ tn'^/N^, which is rapidly decreases with increasing A^. Comparing this amplitude with 
the amplitude of Egrv we get the criterion for when the single flux quantum or M/N flux 




(7.39) 
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quantum periodicity becomes dominant over the fractional ones; this happens when 



However, to investigate this problem we have used the parameter a = ^ <^ 1. The region 
where the single or M/N flux quantum periodicity is dominant is 



One sees that at large L the limits when such behavior exists may be easily recovered. 

In conclusion, then in the dilute density limit or for strong coupling we have discovered 
two different regimes of the Aharonov-Bohm effect, namely, fractional 1/N and integer 
coexisting with fractional M/N and 1/N. Integer and fractional M/N periodic oscillations 
appear simultaneously and may have comparable amplitudes, which are distinct by numerical 
factor. As a rule, in high magnetic fields the amphtude of M/N oscillations is always amaller 
than the amplitude of integer ones. 

Thus, a very surprising implicit role of the Zeeman interaction has been found, i.e. with 
increasing magnetic field, as the contribution of the Zeeman energy increases, in addition to 
the fractional 1/N flux quantum periodicity, there appears the integer- fractional M/N flux 
quantum periodicity. With increasing the magnetic field, i.e. with decreasing M the half 
flux quantum periodicity a,t M — N/2 transforms into fractional M/N ones and integer one. 
Finally, in strong field, when the system is almost polarized, only the integer flux quantum 
periodicity remains. 

This transformation is not continuous but singular. The period of this fine structure 
of the AB effect changes with each spin flip. For almost polarized system, we have the 
coexistence of integer and fractional 1/N flux quantum periodic oscillations, only. On the 




(7.40) 




(7.41) 
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ring with a small number of electrons N < -g^, there occurs the fractional Aharonov- 
Bohm effect. However with increasing A^, when ^ '^^^ single-fractional M/N flux 
quantum periodicity becomes dominant over the fractional 1/A^ effect, although they may 
both coexist. 

VIII. CHANGE OF THE FINE STRUCTURE PERIOD WITH MAGNETIC 

FIELD 

The period of fractional oscillations ^ related to the number of down-spin particles. 
However, with the increase of magnetic field on the ring there may occur spin-flip processes. 
As the result the number of the down-spin particles changes (decreases) and the magneti- 
zation of the ring related to the number M decreases. This Zeeman magnetization may be 
described with the aid of the formula: 

{N - 2M) = [^] (8.42) 

where brackets means an integer part of the value. The value x is a susceptibility of the 
Hubbard chain . First, in thermodynamic limit the value of x has been calculated by 
Takahashi [28] for the case of the half- filling, then for an arbitrary filling the value of x has 



been calculated by Shiba ||29|. For the low density limit a << 1 we discuss the value of x 
may be estimated by the similar way. For the large values of L the thermodynamic limit is 
acceptable. Therefore we use Shiba's result 

whence 
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One sees that with the field the value of M, and, therefore, the period of the fine structure 
decreases 

/r = i-|fl^. (8.45) 

Thus we came to the conclusion that in the Aharonov-Bohm effect of strongly correlated 
system such as the Hubbard ring may occur the oscillation of the smallest amplitude, which 
period decreases, when the field increases. In such an effect the spin orbit interaction might 
be important (see, for example, Ref. [pO|] ) 



IX. TOWARDS REALITY FROM EXACT SOLUTIONS 

In the present work we have presented the exact solution for the Hubbard ring, showing 
a fine structure in the Aharonov-Bohm effect. The mesoscopic rings in existing experiments 
have many channels, large sizes and many electrons ( see, also, discussion in Ref. |^). They 



also have disorder or impurities. The measurements are carried at finite temperatures. Of 
course, all these factors make a large difference between the presented model calculations 
and the realistic situations. 

However on such rings the electron density is much smaller than in a 3D metal, where 
screening exists. Therefore the electron-electron interaction is very important. The effects of 
the strong electron-electron correlations, which give rise the finite structure of AB effect, of 
course, will appear there. But at present time there is no tool to take into account correctly 
the many-body effects related to this interaction. The Bethe ansatz method is only one single 
tool, where the effects of strong correlations may be studied for large number of electrons 
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and a system of a large size, but for some models only. So, using the Bethe ansatz we 
found, that the interaction between particles also may create the finite structure, which is 
characterized by several types of periodic oscillations with the amplitude and period smaller 
than one unit. 

In the integer-fractional regime the effects of disorder and temperature will be to make 
the flux-energy dependence smoother leading to the disappearence of the fine structure. In 
experimental situations, where the disorder or the temperature always exist, but have a 
low level for the ring with large number of electrons one may see three types of oscillations 
in the flux energy dependence. With the next increase of the disorder or the temperature 
the oscillations having the smallest amplitude will disappear. So, for example, instead of 
integer- fractional M/N periodic oscillations on the single ring, one may see only the integer 
flux quantum periodic oscillations. 

The discovered single flux quantum periodicity, when the finite structure has disappeared, 
may be relevant to two recent experimental studies of the Aharonov-Bohm effect on single 
metallic and semiconductor rings Recent experiments [^| have described only the 



full quantum period in the case of individual metallic rings. On a semiconductor single loop 
in the GaAs/GaAlAs system such single flux periodicity has been detected. To explain 
this integer flux quantum periodicity we assume that the individual rings are sufficiently 



disordered. Recently, similar results |]T3[ |^ were obtained that free electrons on the single 
ring with disorder or with the temperature may show only the full single flux quantum 
periodicity. Both these findings are consistent with existing experiments. 

The interaction between particles may also support the full single quantum periodicity, 
and, therefore, it may also be responsible for the observed single flux quantum period. The 
prediction of fractional and fractional M/N oscillations could be made for the ring 
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having lowest level of disorder. Once the resolution of the experimental studies has been 
improved such periodicities may be observed. 

The other experiment which is worth to mention is that the half-flux quantum 

periodicity simultaneously of 10'' rings has been measured. It is also consistent with our 
results. It is reasonable to assume that there is a disorder. So the single ring must show the 
single flux periodicity. For such large ensemble of the rings we must take also into account 
the parity effect. Then it is also reasonable that half of these rings has an even number of 
particles and the other half of these rings has an odd number of particles. Because of the 
discovered and described parity effect, averaging over the two cases of odd and even number 
of particles indicated above, immediately gives the $o/2 period. This interpretation is then 
natural for the experimental results concerning many rings [Q. Thus in the framework of 
our approach we have obtained the unified pictures why the $0 period is seen in individual 
rings and why $o/2 period is seen an ensemble of 10^ rings. 

In addition to these explanations, we predict the existence of oscillations with periods 
smaller than $0 or $o/2- To check the predictions of our theory, it is important to perform 
the Fourier analysis of the persistent current oscillations in order to detect the fractional 
Aharonov-Bohm effect, which coexists with integer one, i.e. whether there are fractional 
oscillations on a background of integer ones. 



Recently the fractional 1 / 4— AB effect has been observed in Auln-iings prepeared by 
e— beam lithography . Since this effect has been observed in the region of a superconducting 
phase transition, where the phase separation occurs. It seems that there the effective system 
reminding the quantum dot chain ring is created. Therefore, the discussed Hubbard model 
for the description of unpaired electrons in these droplets may be applicable. 

In this case the fractional 1 / 4— effect may be created as a result of averaging over rings 
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having different parity of electron's number. This occurs only at a small magnetic field 
when Zeeman energy may be neglected. With the increasing magnetic field the electrons 
become partially polarized and such periodicity is destroyed. But here the fractional M/N 
may appear instead. 

Thus, our finding of a coexistence of the fractional 1/N effect with fractional 1/4 or 1/2 
integer AB effects may not only to explain the observed 1/4 flux quantum periodicity but 
also predicts smallest, namely 1/N and M/N fractional AB effects, once the experimental 
resolution has been improved. 
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Figure Captions 



Fig.l The behavior of the ground state energy as a function of flux / for 6 electrons 
with the single up spin at the values L = 1000 and U = 10 in the region of flux within 
the single fundamental flux quantum. The energy is expressed in the units tlO^ . The zero 
energy corresponds to — 11.9993^. 

Fig. 2 The behavior of the ground state energy as a function of flux / for 11 electrons 
at the values L = 1000 and t/ = 10 in the region of flux within the single fundamental 
flux quantum. One electron has down spin i.e. M — 1. The zero energy corresponds to 
-21.995t 

Fig. 3 The behavior of the ground state energy as a function of flux / for 25 electrons 
at the values L = 7500 and [/ = 10 in the region of flux within the single fundamental flux 
quantum. 12 particles have up-spin and 13 particles have down-spin. The energy is expressed 
in the units tlO^ . The zero energy corresponds to — 24.999544^. 

Fig. 4 The behavior of the ground state energy as a function of flux for 10 electrons 
at the values L = 1000 and [/ = 10 in the region of flux within the single fundamental 
flux quantum. Five particles have up-spin and five particles have down-spin. The energy is 
expressed in the units tlO^ . The zero energy corresponds to — 9.9983t. 

Fig. 5 The behavior of the ground state energy as a function of flux for the ring with 
even number of electrons A'^ and with odd number for projection of angular momentum M, 
for the parameters L = 10^ and U = 10. We present the cases when and M are equal to: 
a) N — 20 and M — 3 , energy is expressed in arbitrary units, the region of flux is the half 
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of fundamental flux quantum, h) N — 30 and M — 5 , energy in units tlO^, the zero energy 
corresponds to — 59.99911t; the region of flux is the half of fundamental flux quantum, c) 
N — 4:0 and M = 7 , energy is expressed in the units tlO^, the zero energy corresponds to 
— 79.997896t; the region of flux is the half of fundamental flux quantum. 

Fig. 6 The behavior of the ground state energy as a function of flux for the ring with 
odd number of electrons N and with even number for projection of angular momentum M, 
for the parameters L — and U — 10. We present the cases when N and M are equal to: 
a) N — 19 and M = 6 , energy is expressed in the units tlO^, the zero energy corresponds 
to — 37.97752t; the region of flux is the half of fundamental flux quantum, b) = 49 and 
M — 6 , energy is expressed in the units tlO"^, the zero energy corresponds to — 97.996131t; 
the region of flux is the half of fundamental flux quantum, c) N — 99 and M — 6 , energy is 
expressed in the units tlO^, the zero energy corresponds to — 197.96808t; the region of flux 
is the quarter of fundamental flux quantum, d) = 99 and M — 23 , energy is expressed 
in the units tlO^, the zero energy corresponds to — 197.96809^; the region of flux is the half 
of fundamental flux quantum. 

Fig. 7 The behavior of the ground state energy as a function of flux for the ring with 
even number of electrons and with even number of a projection of the angular momentum 
M, for the parameters L = 10^ and U = 10. We present the case when A*" = 20 and M = 6, 
energy is expressed in the units tlO^, the zero energy corresponds to — 39.97377t; the region 
of flux is the half of fundamental flux quantum. 

Fig. 8 The behavior of the ground state energy as a function of flux for the ring with 
odd number of electrons N and projection of angular momentum M, for the parameters 
L = 10^ and U — 10. We present the case when A^ = 19 and M — 5, energy is expressed in 
the units tlO^, the zero energy corresponds to — 39.97752t; the region of flux is the half of 
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fundamental flux quantum. 



